Compensation of wake-field-driven energy spread in Energy Recovery Linacs 
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Energy Recovery Linacs provide high-energy beams, but decelerate those beams before dumping 
them, so that their energy is available for the acceleration of new particles. During this deceler- 
ation, any relative energy spread that is created at high energy is amplified by the ratio between 
high energy and dump energy. Therefore, Energy Recovery Linacs are sensitive to energy spread 
acquired at high energy, e.g. from wake fields. One can compensate the time-correlated energy 
spread due to wakes via energy-dependent time-of-flight terms in appropriate sections of an Energy 
Recovery Linac, and via high-frequency cavities. We show that nonlinear time-of-fiight terms can 
only eliminate odd orders in the correlation between time and energy, if these terms are created by 
a beam transport within the linac that is common for accelerating and decelerating beams. If these 
two beams are separated, so that different beam transport sections can be used to produce time-of- 
flight terms suitable for each, also even-order terms in the energy spread can be eliminated. As an 
example, we investigate the potential of using this method for the Cornell x-ray Energy Recovery 
Linac. Via quadratic time-of-flight terms, the energy spread can be reduced by 66%. Alternatively, 
since the energy spread from the dominantly resistive wake fields of the analysed accelerator is ap- 
proximately harmonic in time, a high-frequency cavity could diminish the energy spread by 81%. 
This approach would require bunch-lengthening and recompression in separate sections for acceler- 
ating and decelerating beams. Such sections have therefore been included in Cornell's x-ray Energy 
Recovery Linac design. 



I. INTRODUCTION 

Energy Recovery Linacs (ERLs) accelerate high- 
current particle beams to high energy in a linac. These 
are then used in x-ray l2| , FEL [3j or nuclear physics 
[HHI experiments. Subsequently, the beams are sent into 
the same linac at a decelerating phase to recover the par- 
ticles' energy. This energy is then used to accelerate new 
bunches of particles Q . Only with such energy-recycling 
does it become feasible to accelerate high-current beams 
to high energies in a linac. Today's high-current, high- 
energy beams are produced by storage rings, which reuse 
energetic electrons for millions of turns. The beam emit- 
tance is then limited by the equilibrium emittance that 
is established during these turns. ERLs use each bunch 
of electrons only once, and therefore have the potential 
of producing significantly smaller emittances. 

One of the problems ERLs face is the wake-field-driven 
energy spread that builds up during a pass through the 
ERL. The energy spread at high velocity is relatively 
small, but relevant because it limits the bandwidth of the 
x-ray radiation. More importantly, it is multiplied during 
deceleration by the ratio between high energy and dump 
energy, which is approximately 500 in the case of the Cor- 
nell ERL. Decelerated bunches must be decoupled from 
accelerated bunches, via magnetic fields, in a demerger 
region before the dump. Here, particles with too large 
energy are not bent sufficiently to reach the dump; sim- 
ilarly, particles with too little energy hit the beam pipe 
before the dump, or are even decelerated to zero energy 
before leaving the superconducting linac, and are then 
lost in the cryogenic environment. Therefore, the energy 
spread that bunches have at the end of their deceleration 
has to be limited. 

To illustrate the proposed techniques of wake-field 




FIG. 1: Cornell ERL layout. 1) Injector, 2) linac A, 3) turn- 
around loop, 4) linac B, 5) south x-ray beamlines, 6) CESR, 
7) north x-ray beamlines, 8) first beam dump, 9) second beam 
dump and 10) distributed cryoplant. 



compensation, we use the Cornell ERL, shown in Fig. [T] 
Bunches enter the ERL via the injector at (1). They are 
accelerated first by linac A at (2) and then by linac B 
at (4) after traversing the turn-around loop connecting 
the linacs at (3). They then traverse the Cornell Elec- 
tron Storage Ring (CESR) at (6), producing x-rays in 
beamline sections (5) and (7). Returning to linac A, the 
bunches are decelerated, passed through the turn-around 
loop, and further decelerated by linac B. Finally, the de- 
merger leads decelerated bunches into the beam dump at 
(8). 

One can describe the propagation of a bunch's lon- 
gitudinal phase-space distribution through the ERL via 
transformations for separate ERL components. In the 
Cornell ERL, shown schematically in Fig. [2j the bunch 
center is to be the design particle with t = and energy 
En. The bunch first traverses linac A, which has fre- 
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Now, the bunch returns to linac A, the sec- 
ond turn-around loop, and linac B, which apply 
Eqs. O, © and © to yield 
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{ t7 ,E r } = {t 6 ,E 6 -AE A COS{uJt ^f A) } 
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where <f) A = a + u5t C E, 



(7) 



FIG. 2: Schematic of the Cornell ERL. 



quency to, mapping the initial longitudinal phase-space 
position {to,Eo} for each particle to 



{< 8 , E s } = {t 7 + T56 T a 2 (£ 7 - E c 7 ) 

+ T5m TA 2(E 7 -E 7 -) 2 1 E r } , 



(8) 



{t 1 ,E 1 } = {t Q ,* + ±E A ™^ + ** ) }. (1) 

COS (0a) 

Here, the time coordinate after linac A is chosen so that 
the bunch center arrives at t\ = 0; similarly, we choose 
ti = at the bunch center for each other section. AEa is 
the energy the linac adds to the bunch center. To limit 
the peak accelerating voltage ec ^0 A ) , 4>a will be kept 
below 15°. Next, the particles traverse the first turn- 
around loop with a time-of-flight mapping 



{t 9 ,E 9 } = {t 8 ,E s -AE B C -^±M } , (9) 

COS (0b) 

with 4>' B — 0b +w(5tcE + ^TA2)- ln this state, the bunch 
leaves the ERL for the dump. While energy recovery of 
AEa hi linac A demands that the bunch center return 
to linac A after an odd multiple of half the RF period, 
a possible deviation StcE is included above. Similarly, 
the decelerated beam might require a slightly different 
time than the accelerated beam to pass from linac A to 
linac B, which is described by <5^tA2- 



{ta, E 2 } = {h + T56tai(£i - Ef) (2) 
+ T566 T ai(£i ~E c 1 ) 2 ,E 1 } , 

where the energy of the bunch center is E\ = Eq + AEa- 
Subsequently, linac B applies another accelerating field 
of frequency to to the bunch, mapping {t%, E 2 } to 



II. WAKE FIELDS IN THE CORNELL ERL 

Table U lists the sources and magnitudes of the wake- 
induced energy spread for ERL components. The struc- 
ture of W{t) and its effect on the bunch at CESR for the 
Cornell ERL are shown in Figs. [3] and [4] respectively. 



{t 3 ,E 3 } = {t 2 ,E 2+ AE B COS ^+^ }- (3) 

COS (0b) 

Also 0b is limited to 15°. Next, we add half the accumu- 
lated effect of all wake fields to E 3 : 



{ti, Ei} = {ta, E 3 



W(t 3 ) 
o J ■ 



(4) 



The particles then traverse CESR, with a time of flight 
mapping {£4, E4} to 



{t B , E 5 } = {t A + T56 C e(£ 4 - El) 
+ T566ce(E4-E%) 2 ,E a } 



(5) 



Next, the second half of the accumulated effect of wake 
fields is added: 



{t 6 ,E e } = {t 5 ,E 5 + ^M} . 



(6) 



Source 


Number 


Max (kV/pC) 


7 Cell RF Cavity 


800 


-11.32 


Higher Mode Load (78 mm) 


400 


-0.89 


Higher Mode Load (106 mm) 


400 


-0.50 


Expansion Joint 


356 


-0.74 


Beam Position Monitor (Button) 


664 


-0.35 


Beam Position Monitor (Stripline) 


20 


-0.01 


Flange Joint 


356 


-0.90 


Clearing Electrode 


150 


-0.18 


Gate Valve 


68 


-0.71 


1 m Resistive Wall (12.7 mm) 


2500 


-4.00 


1 m Roughness (12.7 mm) 


2500 


-14.00 


Undulator Taper (3 mm) 


18 


-0.61 


1 m Resistive Wall (3 mm) 


144 


-0.98 


1 m Roughness (3 mm) 


144 


-3.60 



TABLE I: Sources and magnitudes of wake-driven energy 
spread [7|. Listed magnitudes are to be multiplied by the 
bunch charge of 77pC. 



-0,5 



> 
hj 
S 




-2,5 



FIG. 3: Accumulated longitudinal wake potential of the Cor- 
nell ERL. 
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FIG. 5: Illustration of curvature change in E(t) by off-crest 
acceleration and time-of-fiight terms. 



III. COMPENSATION METHODS FOR WAKE 
DRIVEN ENERGY SPREAD 
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FIG. 4: Wake-induced bunch profiles at CESR. Black-top: 
Cosine-like correlated longitudinal phase space from acceler- 
ating on crest with a at = 2ps bunch length. Blue-bottom: 
Longitudinal profile after suffering half the Cornell ERL's 
wake field, ^ . 



A. Time of flight for wake correction 

A bunch's correlated energy spread can be reduced by 
decreasing its slope and curvature in time-energy phase 
space. This can be done by accelerating the bunch off- 
crest in the linacs, choosing an energy-dependent time of 
flight, and decelerating the bunch off-crest. We sketch 
this procedure in Fig. [5J There, arrows show the phase- 
space motion of a particle. Curve 1 shows the initial 
bunch profile; curve 2 shows the bunch after accelera- 
tion by the linacs; curve 3 shows the bunch after applica- 
tion of time-of-flight terms; curve 4 shows the decelerated 
bunch, where the bunch's initial curvature, i.e. second- 
order correlation between time and energy, is eliminated. 
The arrows follow one particle of the bunch through this 
process. 



1. Linear time of flight for linear-wake correction 



We first eliminate the average slope ^S 2 - = 
£9(2) ^3 ( — g_) £ phase-space bunch profile before the 
dump, where T is chosen to be six times as large as crt, 
the rms temporal bunch length. Eliminating A_Eg may 
not necessarily decrease the energy spread, but serves to 
symmetrize the wake for higher-order compensation. We 
minimize |A£Jg| using time-of-flight terms. Additionally, 
T56tai can be used to minimize the energy spread |Ai?4| 
in CESR. This minimization has to be done numerically, 
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and we start with the T56tA2 that eliminates the phase- 
space slope ^ a at the bunch center before the dump. 



We obtain from 

dtg 



dEj 
dti 



dEj 
dtp 
dti 
ill, i 



with and 4^ horn 

dto dto 



(lz; 




dt a 



M 



./0 



dzo 
dt 







(10) 



(11) 



2. Nonlinear time of flight for nonlinear-wake correction 

With the average phase-space slope of the bunch pro- 
file before the dump eliminated, the next step is to re- 
duce the phase-space curvature at the bunch center 
before the dump. Unfortunately, minimizing the abso- 
lute curvature at the bunch center does not necessarily 
minimize the total energy spread. Therefore, we again 
use numerical minimization, which we start with analyt- 
ically determined parameters that eliminate d d ^ 9 with 



J = 



1 

-1 



where we have chosen — for phase-space distri- 
butions that enter the linac without linear time-energy 
correlation. Here, 



j'+i 

M ij = l[M kk _ 1 

k—i 



M 



kk-l 



d{t kl E k ) 
d(tk-i,Ek- 



(12) 



describes the transfer matrices in longitudinal phase 
space; 



M 



1 
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Mai 



-LuAEAtan(4>A) 1 



1 T56 T ai 







1 



(13) 



(14) 



2 d A Ei dtj _ dEj drtt 

d Ei dtf t dt dt dtp 

dtf ~ rftp 
dt 



(21) 



dtj \ -^dz- dfzj _ f dtj \ ~ 3 dzf T j 
dt J dt ~dtl ~\dt Q ) dt ~ l0 ~ 



, d z% 
dti 



The following iteration determines 



d 2 z % 
dti 



M 



d 2 z l - 1 



9i : 
H, 



dzf-i 



H 



dzi-i 



dti 



with 



dto 1 dto 

( I[M«.-i^.-i)l 

[HlEiiU-uEi^!)] 



(22) 
(23) 
(24) 



where H[f(x, y)] is the Hessian matrix of f(x,y); 





ujAE b 




Mo 2 =l 1 ° I , (15) ri[f(x,y)]=[ % % a f I , (25) 

— 32 ^ -wA£ B tan(<£ B ) l) 1 * ' ' ' ' 

dzi 

cLIlU U1C 
r\a +r\ 

M4 3 =M 65 = ( wuo) , ) . (16) 



and the ^ are determined by Eq. (flT]) . Equation ([22]) 
leads to 



^$=yM ik 9k+M i0 [ / Eo I , (2(i! 




with — ttt being the curvature in longitudinal phase space 
with which the bunches enter the linac. We here as- 

(1 \ sume that the bunch enters the ERL without time-energy 

ujAEa sin( -f A \ 1 J ' ^ 18 ' correlation; i.e. ^ = ^Mp- = 0. With Eq. Cm 

"cos(0a) / at o dt 1 ' ' 

and using the symplecticity of transport matrices, i.e. 



M T J M = J, Eq. (26]) leads to 



Msr = : °: TA2 , (19) 



*r = U) (M)Em^- (27) 
/ i o \ 

—98 = ( a sin(0{,) -, ) ■ (20) The Hessian matrices evaluated at the bunch center are 



cos(<^b) / zero except for 
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H[Ei(t ,Eo)] = 



H[t 2 (ti,£?i)] = 
H[E 3 (t 2 ,E 2 )} = 




H[E 4 (t 3 ,E 3 ))=H[E 6 (t 5 ,E 5 )] 




FIG. 6: Illustration of curvature cancellation in E(t) by time- 
(31) of-flight terms in a common turn-around loop for accelerating 
and decelerating beams. The ordinate specifies the energy 
relative to the bunch center in MeV. 



H[t 5 {U,E 4 )] 





2T566ce 



H[E 7 (t e ,E 6 )] 



lu 2 AE a 



H[t 8 {t 7 ,E 7 )} = 



cos(0 A ) 






2T566 T A2 



H[E 9 (t s ,E s )} = 



lj 2 AE b co 



(0b) 



(32) 



(34) 



(35) 



3. CESR time of flight 

As explained in Fig. [5J a first-order energy correlation 
has to be admitted in CESR to influence the curvature 
of the longitudinal distribution. Consequently, the bunch 
entering CESR has a large energy spread, which undesir- 
ably broadens the bunch's x-ray spectrum, rendering it 
unfeasible to use time-of-flight terms in CESR for wake 
compensation. 



4- Common turn-around loop for two ERL beams 

In this section, we choose <j/ A — <fi A and <fi' B — <fi B 
for balanced acceleration and deceleration in each linac. 
Then, for the reference energy of the turn-around loop 
to be the same for both accelerating and decelerating 
beams, E\ = E%, one has to choose AE A = AE B +W(0), 
or for a total acceleration of AE 
AE+W ^ and AE B - AE ~ W ^ 



AE A + AE B , AEp 



The slope, at the bunch center in CESR is to be 
eliminated to have small energy spread for x-ray experi- 
ments. Consequently, = and 



(33) dz 7 



, = M 7 4 dto 
dt — 74 I 



dh / 1 

dt \ ujAEAtaxiU 



W'(0) 



(36) 



so Eq. ([TO]) gives the slope at the bunch center in the 
second pass through the turn-around loop as ^r 1 = 

u)AE A ta,n(<j) A ) + KP-. Eq. ((TUJj also gives the slope at 
the bunch center in the first pass through the turn-around 
loop, = — ujAE a tan(0A), and the beam therefore 
enters the turn-around loop with roughly opposite slopes 
on each pass; ff^ -ff • 

As a result, the net curvature change in the turnaround 
is zero, as explained in Fig. [6l First, linac A accelerates 
off-crest to produce a correlated energy spread here with 
positive slope. Second, the turn-around loop adds a sec- 
ond order time-of-flight shift, here to the right and thus 
effectively to the bottom of the bunch. In the second 
pass through the turn-around loop, the slope is negative, 
so that the second order time-of-flight terms shift effec- 
tively to the top of the bunch, compensating the curva- 
ture change in the first pass. The curvature from the 
first turn-around loop is counterbalanced by the curva- 
ture from the second turn-around loop. Evidently, this 
method cannot change even-order time-energy correla- 
tions. Only odd-order time-energy correlations at the 
dump can be eliminated when is to be eliminated. 



5. Separate turn-around loops for each ERL beam. 

In the previous two sections, we have seen how longitu- 
dinal phase-space curvature reduction can diminish the 
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bunch energy spread, but can be feasibly carried out in 
neither CESR nor a single turn-around loop. However, 
passing the bunch through different turn-around loops 
before and after CESR enables such curvature reduction. 

In this section, we choose <f>' A — 4>a and <fi' B — (f>B 
for balanced acceleration and deceleration in each linac. 
The variable T56tai is chosen to eliminate the slope, ^J^- 
from Eq. (fit)]) , at the bunch center in CESR; 



T56tai = [wAEa tan(</>A 
t W(0) 



>AE B tan((/) B )}- 1 , 
and T56ta2 is chosen to eliminate at the dump 



T56 TA 2 = -[^ASBtan(^B)]" 1 



(38) 



r W'(0) 



;A£ A tan^)]" 1 



5000 



4999.5 



(37) S 




Similarly, T566tai will be used to eliminate the curva- 
ture, from Eq. (J2TJ) , at the bunch center in CESR, 
and T566tA2 to eliminate that at the dump, leading to 



T566tai 



T566ta2 



2W"(0) -Alo 2 AE b 
[2LuAE B tan((j) B ) - W'iO)] 3 
AE A 

2aj[AE A t&n((f)A)} 3 ' 



2W"(0) + Au 2 AE A 
[2wA£ A tan(0A) + W'{0)} 3 
AE B 



(39) 



(40) 



2w[A£ B tan(0 B )] 3 



This choice compensates the curvature of the wake, as 
well as the curvature of the RF cosine wave. For the 
Cornell ERL, W'(Q) ~ 0, simplifying Eqs. ((Ml) and (|30|). 
The part in Eqs. (f3"9")l and (|4l7|) that remains when setting 
^'(0) = W"(0) = corrects the curvature of the accel- 
eration and deceleration cosine function Q ; the rest is re- 
sponsible for correcting the curvature of the wake fields. 
The phases 4>a and 0b can be chosen independently Sub- 
sequently we choose 0b so that T56tai = T^Q^A2 = 
when no wake field is present; i.e. AEa tan(0A) = 
A£ B tan((/>B). 

We then numerically find the combination of T56tai 
and T566tai that minimizes the maximum energy dif- 
ference between bunch particles at CESR, starting with 
Eqs. (|3T|) and (f3"9"| . Subsequently, we numerically find the 
combination of T56tA2 and T566tA2 that minimizes the 
maximum energy difference between bunch particles at 
the dump, starting with Eqs. (|38[) and ([40]) . As an exam- 
ple we use AE A = AE B = 2495MeV and cj) A = -0b = 
15°. Figure [7] shows the best result, which has been ob- 
tained for T56 T ai = 4.4 x lO^psMeV- 1 , T566 T ai = 



FIG. 7: Bunch profiles after nonlinear-wake correction with 
separate turn-around loops for each ERL beam. Left ordi- 
nate gives energy at CESR while right ordinate gives energy 
at dump. Black-top: Cosine-like correlated longitudinal phase 
space from accelerating on crest with a at = 2ps bunch length. 
Blue-middle: Longitudinal profile after suffering half the Cor- 
nell ERL's wake field. Red-bottom: Longitudinal profile at 
dump with residual energy spread due to higher-order corre- 
lations. 



7.7 x lCT 4 psMeV- 2 , T56 T a 2 = -5.6 x lO^psMeV" 1 , 
T566 T A2 = 1-5 x lCT 3 psMeV- 2 and T56 C e = T566 C e = 
0, with the maximum energy difference between bunch 
particles decreasing by 66%, from 3.2 to l.lMev. The 
residual energy spread is due to higher-order terms in 
W(t). Hence, we acquire a large decrease in energy 
spread, without having had to send the bunch into CESR 
with an undesirable energy spread. 



B. Harmonic- wake correction 

1. High-frequency cavity 

Instead of attempting to minimize phase-space curva- 
ture at the bunch center before the dump by nonlinear 
time-of-fiight terms, one could add energy to the bunch 
to flatten its longitudinal phase-space distribution. Wc 
simulate this by inserting a cavity with period T ap- 
proximately six times as large as crj, the rms temporal 
bunch length, immediately after CESR. T is chosen to 
ensure the cavity is a multiple of the linac frequency, so 
subsequent bunches have the same phase at the cavity. 
The cavity maps the phase-space position after CESR, 
{t e ,E 6 } from Eq. ©, to 



7 



{t' 6 ,E' 6 } = {t 6 ,E 6 + AE 1 cos( 



2n[t 6 



T 



(41) 



5000 



so *6 and E 6 in Eq. ([7|) become t' 6 and respectively. 

In Eq. (1411) . t* and AE± are determined from Fourier 
analysing Wif): 



oo 



07777- 

lF(t)-^A£;„cos(^[t 5 

n=0 



CD 



(42) 



We again choose A = 0a and 4>' B = 0b for balanced 
acceleration and deceleration in each linac, and we choose 
0b to satisfy Eq. ([57]) . which eliminates the slope at the 
bunch center in CESR. 

Using a single turn-around loop, we numerically find 
the combination of T56ta = Tbd^Ai = T56^ A2 , t* 
and AEi that minimizes the maximum energy differ- 
ence between bunch particles at the dump, starting with 
Eq. (J42|) . As an example we use AEa — AEb — 
2495MeV and a = 10°, which is non-zero because 
there has to be a slope in longitudinal phase space to 
eliminate the first-order correlation via T56ta- Fig- 
ure [5] shows the best result, which has been obtained 
for t\ = -9.0 x 10- 2 ps, AE 1 = 1.5MeV, T56 TA = 
-1.3 x lO^psMeV" 1 , B = -9.9°, and T566 T ai = 
T566tA2 = T56ce — T566ce = 0, with the maximum 
energy difference between bunch particles decreasing by 
81%, from 3.2 to 0.60MeV. The residual energy spread is 
due to higher harmonics in W(t). The decrease in energy 
spread is larger than that obtained using nonlinear-wake 
correction, indicating the effectiveness of harmonic cor- 
rection. 

Although the decrease in energy spread is substantial, 
the required cavity frequency l£ ~ 80GHz is too high 
to be economically feasible. The energy that has to be 
added to the bunch center is approximately 2MV, for 
a beam current of 0.1 A. While power sources of up to 
40GHz could be made available, a power of 200kW seems 
unfeasible. An X-band power sources at approximately 
11.7GHz might be feasible, but the bunch length would 
have to be increased to do harmonic- wake correction with 
this lower frequency. 



2. Bunch expansion and recompression 

a. Common turn-around loop for two ERL beams. 
Reducing the frequency to 11.7GHz requires expanding 
the 6<7t of the bunch by 7 from 12ps to 84ps. Because a 
longer bunch acquires more energy spread in the cosine- 
like RF field, we should recompress the bunch immedi- 
ately after the 11.7GHz cavity. 

Resizing the bunch in CESR requires energy spread 
in CESR, and is therefore unfeasible. Using the turn- 
around loop allows for the compensation of energy spread 
due to 4>a via a suitable 0b • Here we simulate a cavity in 
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FIG. 8: Bunch profiles after harmonic- wake correction with 
high-frequency cavity. Left ordinate gives energy at CESR 
while right ordinate gives energy at dump. Black-top: Cosine- 
like correlated longitudinal phase space from accelerating on 
crest with a a t = 2ps bunch length. Blue-bottom: Longitu- 
dinal profile after suffering half the Cornell ERL's wake field. 
Red-middle: Longitudinal profile at dump with residual en- 
ergy spread due to higher harmonics. 



the center of a single turn-around loop, with T56ta de- 
scribing the first half of the turn-around loop, to map 
the phase-space positions, {£2, £-2} and {tg.Eg,} from 
Eqs. © and (©, to 

{4 E' 2 } = {t 2 , E 2 - AE X cos( 27r[h ~ tt] )} , (43) 
{4 E>} = {t 8 , E s + AE, cos( 27r[f8 r -% , (44) 

where t\ and AE\ are determined from Eq. p2|) . with 
W(t) replaced by the difference between £7 and Ef\ 
which is E7 obtained from setting W(t) = 0. Here, T 
is approximately six times the new rms temporal bunch 
length. T is chosen to ensure the cavity frequency is an 
odd multiple of the linac frequency, so the bunch has the 
opposite phase on both passes through the cavity. The 
second half of the turn-around loop maps {t^E^} and 
{f 8 ,££} to 



{t'J,.E'l} = {t' 2 + TW TA {E' 2 -E' 2 c ) 

+ T566' TA (E 2 -E 2 C ) 2 ,E> 2 } , 

{tlE'i} = {t' s + Th&' TA (E' % ~E'£) 

+ T5m' TA (E' 8 -E> 8 c f,E / 8 } , 



(45) 
(46) 
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where we choose T56' TA — — T56ta to correct the bunch 
length change from the first half of the turn-around 
loop. The numerical minimization will lead to a devi- 
ation of these two time-of-flight terms in order to find 
the smallest energy spread. For simplicity, we assume 
T566^ A = — T566ta = 0, and do not use time-of-flight 
terms in CESR. Eqs. (gSJ) and (gBJ) imply that t, and Ei 
in Eqs. and (0 become t'( and E" respectively, for 
i = 2,8. 

We again choose <j)' A = 0a and (f)' B = 0b for balanced 
acceleration and deceleration in each linac, and 0b to 
eliminate the slope at the bunch center in CESR. 

To determine the T56ta suitable for this method, we 

dig 

obtain the expansion factor = |§- , using and 4J 1 
from Eq. Then, at the bunch center, 



= 1 - T56tawAS b tan(0 B ) , (47) 

dt-j 

where W(t) and AE± have been set to zero. To ex- 
pand the decelerating beam by a factor of 7, we require 

T56 T A = - w AB B tan(0 B )- 

During the first pass through the turn-around loop, we 
obtain 



-1 = 1- T56 T aojAE a tan(0 A ) = -5 , (48) 
at i 

with Ai?A tan(0A) = — AEb tan(0 B ). The bunch is 
therefore expanded in the high-frequency cavity during 
both passes through the turn-around loop. 

We then numerically find the combination of T56ta', 
t\ and AE\ that minimizes the maximum energy dif- 
ference between bunch particles at the dump, starting 
with Eq. (H2). As an example we use AEa — A_E B — 
2495MeV, A = 15°. Figure M shows the best result, 
which was obtained for t\ = 5.9 x 10~ 2 ps, A£?i = 
0.14MeV, T56 T a' = -LlpsMeV" 1 and B = -15°, with 
the maximum energy difference between bunch particles 
decreasing by 78%, from 3.2 to 0.72MeV. 

This result is decent, but the method is unfeasible be- 
cause the bunch is over-compressed during the first pass, 
and therefore has to pass a point where it is fully com- 
pressed to very short length. At this point, Coherent 
Synchotron Radiation (CSR) damage to the beam can 
be very severe. 

b. Separate turn-around loops for each ERL beam. 
The problem of over-compression and therefore very 
short bunches with strong CSR effect can be overcome by 
passing the bunches through different turn-around loops 
before and after CESR, such that only one of the turn- 
around loops carries out harmonic-wake correction. The 
harmonic correction is performed in the turn-around loop 
carrying the bunch after CESR, while the turn-around 
loop carrying the bunch before CESR has T56tai chosen 
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FIG. 9: Harmonic- wake correction with bunch expansion and 
recompression, in a single turn-around loop. Left ordinate 
gives energy at CESR while right ordinate gives energy at 
dump. Black-top: Cosine-like correlated longitudinal phase 
space from accelerating on crest with a at = 2ps bunch length. 
Blue-bottom: Longitudinal profile after suffering half the Cor- 
nell ERL's wake field. Red-middle: Longitudinal profile at 
dump. 

according to Eq. (j3"T|) to eliminate the slope at the bunch 
center in CESR. The same expansion factor, = 7, and 
parameter choices, with the exceptions of t\ = — 0.61ps, 
AEi — 0.21MeV, allow us to decrease the maximum en- 
ergy difference between bunch particles by 78%, from 3.2 
to 0.69Mev. The result is shown in Fig. ITD1 

IV. CONCLUSION 

We have investigated the potential of using time-of- 
flight terms in ERL loops to reduce wake-driven corre- 
lated energy spread. As an example, we have used the 
Cornell ERL. One could use the time of flight in the 
ERL return pass that contains the x-ray sources, but this 
method requires undesirable correlated energy spread in 
this pass. We found that time-of-flight terms in the linac 
that are common to the accelerating and the decelerating 
beams cannot compensate even orders in the time-energy 
correlation. However, when the accelerating and the de- 
celerating beams are separated in the linac, time of flight 
for the two passes through the linac can be used to re- 
duce the energy spread in the return loop, as well as in 
the beam dump. The correlated energy spread is reduced 
by 66% in our example. 

Furthermore, we have looked into the potential of us- 
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FIG. 10: Harmonic-wake correction with bunch expansion 
and recompression, in separate turn-around loops for each 
ERL beam. Left ordinate gives energy at CESR while right 
ordinate gives energy at dump. Black-top: Cosine-like corre- 
lated longitudinal phase space from accelerating on crest with 
a at = 2ps bunch length. Blue-middle: Longitudinal profile 
after suffering half the Cornell ERL's wake field. Red-bottom: 
Longitudinal profile at dump. Bunch profiles are similar to 
those in Fig. [8] for harmonic-wake correction without bunch 
expansion and recompression. 



ing high-frequency cavities to reduce the energy spread. 
We find that placing an 80Ghz cavity after CESR could 
eliminate the first harmonic of the energy spread, and de- 
creases the energy variation by 81% in our example. Un- 
fortunately, a cavity with such a high frequency is unfea- 
sible, necessitating bunch-lengthening to permit a lower 
cavity frequency. The method we have analysed involves 
expanding the bunch in the first half of the turn-around 
loop, correcting the first harmonic of the time-expanded 
energy spread in a cavity, then recompressing the bunch 
in the second half of the turn-around loop. Separate turn- 
around loops would be needed to avoid strong CSR forces 
during overcompression, yielding a 78% decrease in en- 



ergy variation in our example. 

Since the wake-induced energy spread is better de- 
scribed by a sinusoid than a parabola, harmonic correc- 
tion is somewhat more effective than nonlinear time-of- 
flight terms at minimizing the energy spread at the dump. 
But bunch compression and decompression, as well as a 
section with high frequency structures, are hard to de- 
vise. In both cases, two turn-around loops are required, 
and have thus been included in the Cornell ERL design. 



> 



9, 5 



9 



3, 5 



7 , 5 







\ 

\ 


/ ; 




y ] % 


I \ \ 


r \ 

f \ 

/ 



□ 2 
t(ps) 



FIG. 11: Results. Black-top: Longitudinal profile at dump 
without wake-correction. Blue-middle: Dump profile with 
harmonic wake-correction. Red-middle: Dump profile with 
nonlinear time-of-flight wake-correction. Harmonic wake- 
correction reduces energy spread more but is less feasible than 
nonlinear wake-correction. 
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